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Abstract. We study thresholds for the appearance of a 2-core in random 
hypergraphs that are a mixture of a constant number of random uniform 
hypergraphs each with a linear number of edges but with different edge 
sizes. For the case of two overlapping hypergraphs we give a solution 
for the optimal (expected) number of edges of each size such that the 
2-core threshold for the resulting mixed hypergraph is maximized. We 
show that for adequate edge sizes this threshold exceeds the max;imum 
2-core threshold for any random uniform hypergraph, which can be used 
to improve the space utilization of several data structures that rely on 
this parameter. 

1 Introduction 

The 2-core of a hypergraph H is the largest induced sub-hypergraph (possibly 
empty), that has minimum degree at least 2. It can be obtained via a simple 
peeling procedure (Algorithm 1) that successively removes nodes of degree 1 
together with their incident edge. 

Algorithm 1: Peeling 
Input: Hypergraph H 

Output: Maximum induced sub-hypergraph with minimum degree 2. 
while H has a node v of degree < 1 do 

if V is incident to an edge e then remove e from H 

remove i; from H 

return H 

Let i?^ p be a random fc-uniform hypergraph with n nodes where each of the 
possible (^') edges is present with probability p independent of the other edges. 
In the case that the expected number of edges equals c • n for some constant 
c > 0, the following theorem (conjectured e.g. in [17], rigorously proved in [19] 
and independently in [13[) gives the threshold for the appearance of a 2-core in 
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where Po[A] denotes a Poisson random variable with mean A. 
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Theorem 1 ([19, Theorem 1.2]). Let k > 3 be constant, and let c*{k) = 
min;^>o t{X, k). Then for p = c ■ n/ with probability 1 — oil) for n — > cxd the 
following holds: 

(i) if c < c* then ^ has an empty 2-core, 
(a) if c > c* then H!^ ^ has a non-empty 2-core. 

Remark 1. Actually this is only a special case of [19, Theorem 1.2] which covers 
£-cores for ft-uniform hypergraphs for all £ > 2, k and £ not both equal to 2. 

Now consider a mixture of graphs ^ on n nodes for different values of p and k. 
Let H^p be a random hypergraph with n nodes where each of the possible (^ ) 
edges is present with probability pi, given via the vectors k — (fci, fc2, . . . , kg) 
and p = {pi,P2t ■ ■ iPs)- While studying cores of hypergraphs in the context 
of cuckoo hashing the authors of [7] described how to extend the analysis of 
uniform hypergraphs to mixed hypergraphs, which directly leads to the following 
theorem. For a. = (ai, q;2, . . . , ctg) G [0, 1]* with oti — ^ l^t 

*(^' «) = TIP! ■ (2) 

E • fc. • ( Pr (Po [A] > 1) 

Theorem 2 (generalization of Theorem 1, implied by [7]). Let s > 1 

he constant. For each 1 < i < s let ki > 3 be constant, and let ai G [0, 1] 
be constant, where X]i=i ~ 1- Furthermore let c*{k,a) = minA>o ^(A, fc, a). 
Then for pi — ai ■ c ■ n/ {J^^ with probability 1 — o(l) for n — > oo the following 
holds: 

(i) if c < c* then p has an empty 2-core, 
{a) if c > c* then p has a non-empty 2-core. 

Using ideas from [7, Section 4] this theorem can be proved along the lines of [19, 
Theorem 1.2[ utilizing that H^p is a mixture of a constant number of indepen- 
dent hypergraphs. 

Remark 2. Analogous to Theorem 1, Theorem 2 can also be generalized such 
that it covers ^-cores for all £ > 2. 

Now consider hypergraphs H'^ p with edge probabilities pj = a, • c • n/(^ ) as in 
Theorem 2. One can ask the following questions. 

1. Assume k is given. What is the optimal vector a* such the that threshold 
c* {k, ex*) —: c* (k) is maximal among all thresholds c* (fe, a)? In other words, 
we want to solve the following optimization problem 

c*(fc) = mini(A, k, a*) = maxmint(A, fc, a) . (3) 

A>0 a A>0 

2. Is there a k such that a* gives some c*(fe) that exceeds the maximum 
2-core threshold c*(fc) among all fc-uniform hypergraphs (not mixed), which 
is known to be about 0.818 for fc = 3, see e.g. [12,17, conjecture[, [5, proof[. 

Remark 3. Often the 2-core threshold is given for hypergraph models slightly 
different from p. The justification that some "common" hypergraph mod- 
els are equivalent in terms of this threshold is given in Section 1.2. 



2 



1.1 Results 

We give the solution for the non-hnear optimization problem (3) for s — 2. That 
is for each k = {ki, we either give optimal solutions a* = {a* , 1 — a*) and 
c* (fc) in analytical form or identify a subset of the interval (0, 1] where we can 
use binary search to determine a* and therefore c*{k) numerically with arbitrary 
precision. Interestingly, it turns out that for adequate edge sizes ki and fc2 the 
maximum 2-core threshold c*{k) exceeds the maximum 2-core threshold c*{k) 
for A:-uniform hypergraphs. The following table lists some values. 



(fcl, ^2) 


(3,3) 


(3,4) (3,6) (3,8) (3,10) (3,12) (3,14) (3,16) 


(3,21) 


c* 
a* 
k 


0.81847 
1 
3 


0.82151 0.83520 0.85138 0.86752 0.88298 0.89761 0.91089 
0.83596 0.85419 0.86512 0.87315 0.87946 0.88464 0.88684 
3.16404 3.43744 3.67439 3.88795 4.08482 4.26898 4.47102 


0.92004 
0.88743 
5.02626 



Table 1. Optimal 2-core thresholds c*(fe), k = (fci,fc2), and a* — (a*,l — a*), and 
A; = a* ■ fci + (1 — a*) • ^2- The values are rounded to the nearest multiple of 10~^. 

More comprehensive tables for parameters 3 < fci < 6 and ki < k2 < 50 are given 
in Appendix B. The maximum threshold found is about 0.92 for k — (3, 21). 

Remark 4- So why does it help to use edges of different sizes? Consider a k- 
uniform hypergraph H^^ ^ that has a non-empty 2-core with node set V and edge 
set E. Let V' be the set of nodes outside V, that is V' D V = ^. Assume that c 
is just above c*(fc). Then there are many small sets Ei,E2, ■ ■ ■ C E, such that if 
one removes all edges of any of these sets, the 2-core of the remaining hypergraph 
would be empty. Now randomly replace a constant fraction /3 = 1 — a of the edges 
of p by edges of larger size. If /3 is large enough, then it is likely that there 
exists a set Ei , where all of the edges are substituted and all of the corresponding 
larger edges are incident with nodes from V' . Consider an arbitrary large edge 
e with e n y 7^ 0. If /? is small enough, then it is likely that there is at least one 
node V from eOV that is incident to only one large edge (namely e) . It follows 
that e will be removed by the standard peeling algorithm (Algorithm 1). Hence, 
if /3 is not too small and not too large, then it is likely that there exists a set Ei 
whose edges are substituted by larger edges that will be removed by the peeling 
algorithm, which results in an empty 2-core. 

1.2 Extensions to Other Hypergraph Models 

While Theorems 1 and 2 are stated for hypergraphs H!^p, one often considers 
slightly different hypergraphs, e.g. in the analysis of data structures. 

Let H^ „i,a ^iid Hn,m,a t>e random hypergraphs with n nodes and m edges, 
where for each 1 < i < s, a fraction of of the edges are fully randomly 
chosen from the set of all possible edges of size ki. In the case of i?,^„j „ the 
random edge choices are made without replacement and in the case of „ ^ 
the random edge choices are made with replacement. Using standard arguments, 
one sees that if m = c ■ n, ki > 3, and pi = (Xi ■ c - nj (^) as in the situation of 
Theorem 2, the 2-core threshold of El^ra a the same as for H^p (see e.g. [9, 
analogous to Proposition 2] ) , and the 2-core threshold of „j ^ is the same as 
for H!^jn,a (see e.g. [9, analogous to Proposition 1]). 
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1.3 Related Work 



Non-uniform hypergraphs have proven very useful in the design of erasure cor- 
recting codes, such as Tornado codes [16,15], LT codes [14], Online codes [18[, and 
Raptor codes [20[ . Each of these codes heavily rely on one or more hypergraphs 
where the hyperedges correspond to variables (input/message symbols) and the 
nodes correspond to constraints on these variables (encoding/check symbols). An 
essential part of the decoding process of an encoded message is the application of 
a procedure that can be interpreted as peeling the hypergraph (see Algorithm 1) 
associated with the recovery process, where it is required that the result is an 
empty 2-core. Given m message symbols, carefully designed non-uniform hyper- 
graphs allow, in contrast to uniform ones, to gain codes where in the example of 
Tornado, Online, and Raptor codes a random set of (1 -|-e) • m encoding symbols 
are necessary to decode the whole message in linear time (with high probability) , 
and in the case of LT codes a random set of m + o(m) encoding symbols are 
necessary to decode the whole message in time proportional to m ■ ln(m) (with 
high probability). Tornado codes use explicit underlying hypergraphs designed 
for a given fixed code rate, whereas LT codes and its improvements. Online and 
Raptor codes, use implicit graph constructions to generate essentially infinite 
hypergraphs resulting in so called rateless codes. In the case of Torndado codes 
the size of the hyperedges as well as the degree of the nodes follow precalculated 
sequences that are optimized to obtain the desired properties. In the case of LT 
codes, as well as in the last stage of Raptor and Online codes each node chooses 
its degree at random according to some fixed distribution, and then selects its 
incident hyperedges uniformly at random. (For Online codes also a skewed selec- 
tion of the hyperedges is discussed, see [18, Section 7[.) While the construction 
of the non-uniform hypergraph used for these codes is not quite the same as for 
m a (or H^p, H!^ ma)' sincc, among other reasons, the degree of the nodes is 
part of the design, they are similar enough to seemingly make the optimization 
methods / heuristics of [15[ applicable, see footnote [11, page 10[. Having said 
that, compared to e.g. [15[, our optimization problem is easier in the sense that 
it has fewer free parameters and harder in the sense that we are seeking a global 
optimum. 

1.4 Overview of the Paper 

In the next section we discuss the effect of our results on three succinct data 
structures. Afterwards, we give our main theorem that shows how to determine 
optimal 2-core thresholds for mixed hypergraphs with two different edge sizes. 
It follows a section with experimental evaluation of the appearance of 2-cores 
for a few selected mixed hypergraphs, which underpins our theoretical results. 
We conclude with a short summary and an open question. 

2 Some Applications to Succinct Data Structures 

Several succinct data structures are closely related to the 2-core threshold of 
fc-uniform hypergraphs (often one considers ^ a ^'^^ s — ai — 1 and fci = 3). 
More precisely, the space usage of these data structures is inversely proportional 
to the 2-core threshold c*(fc), while the evaluation time is proportional to the 
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edge size k. By showing that the value of c* (3) can be improved using mixed 
hypergraphs instead of uniform ones, our result opens a new possibility for a 
space-time tradeoff regarding these data structures, allowing to further reduce 
their space needs at the cost of a constant increase in the evaluation time. Below 
we briefly sketch three data structures and discuss possible improvements, where 
we make use of the following definitions. 

Let a = (ai, a2, . . . , a„) be a vector with n cells each of size r bits. Let 
S = {xi, X2, ■ ■ ■ , Xm} be a set of m keys, where S is subset of some universe U 
and it holds m = c • n for some constant c < 1. The vector cells correspond to 
nodes of a hypergraph and the keys from S are mapped via some function ip to 
a sequence of vector cells and therefore correspond to hyperedges. We identify 
cells (and nodes) via their index i, 1 < i < n, whereas stands for the value of 
cell i. The following three data structures essentially consist of a vector a and a 
mapping tp. For each data structure we compare their performance, depending if 
if realizes a uniform or a mixed hypergraph. In the case of a uniform hypergraph 
each key Xj is mapped to fc = 3 random nodes (p{xj) = {gi(xj), g2{xj)^ g^lxj)) 
via functions gi,g2,93 ■ U —?' {1,2, ... ,n}. In the case of a mixed hypergraph, 
as an example, a fraction of a* — 0.88684 keys are mapped to 3 random nodes 
using functions gi, g2, gs and a fraction of 1 — a* keys are mapped to 16 random 
nodes via functions g'i,g2,. .. ,g'iQ :f7— ?>{l,2,...,n}. We fix c below the 2-core 
threshold to c = c* — 0.005, which gives c — 0.813 in the uniform case and 
c = 0.906 in the mixed case, cf. Table 1. The reason why we use a rather small 
distance of 0.005 is that for large m one observes a fairly sharp phase transition 
from "empty 2-core" to "non-empty 2-core" in experiments, cf. Section 4. 

2.1 Invertible Bloom Lookup Table 

The invertible Bloom Lookup Table [11] (IBLT) is a Bloom filter data structure 
that, amongst others, supports a complete listing of the inserted elements (with 
high probability). We restrict ourselves to the case where the IBLT is optimized 
for the listing operation and we assume without loss of generality that the keys 
from S are integers. Each vector cell contains a summation counter and a quan- 
tity counter, initialized with 0. The keys arrive one by one and are inserted into 
the IBLT. Inserting a key Xj adds its value to the summation counter and in- 
crements the quantity counter at each of the cells given via ip{xj). To list the 
inserted elements of the IBLT one essentially uses the standard peeling process 
for finding the 2-core of the underlying hypergraph (see Algorithm 1). While 
there exists a cell where the quantity counter has value 1, extract the value of 
the summation counter of this cell which gives some element Xj . Determine the 
summation counters and quantity counters associated with Xj via evaluating 
f{xj) and subtract Xj from the summation counters and decrement the quan- 
tity counters. With this method a complete listing of the inserted elements is 
possible if the 2-core of the hypergraph is empty. Therefore in the case of uni- 
form hypergraphs we get a space usage of n/ra ■ r w 1.23 • r bits per key. As 
already pointed out by the authors of [11], who highlight parallels to erasure 
correcting codes (see Section 1.3), a non-uniform version of the IBLT where keys 
have a different number of associated cells could improve the maximum fraction 
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c = m/n where a complete listing is successful with high probability. Using our 
example of mixed hypergraphs leads to such an improved space usage of about 
1.10 • r bits per key. 

2.2 Retrieval Data Structure 

Given a set of key- value pairs {[xj,Vj) \ Xj G S,Vj e R,j € [m]}, the retrieval 
problem is the problem of building a function f : U ^ R such that for all Xj from 
S it holds f{xj) — Vj\ for any y from U \ S the value f{y) can be an arbitrary 
element from R. Chazelle et al. [4] gave a simple and practical construction of 
a retrieval data structure, consisting of a vector a and some mapping that 
has constant evaluation time, via simply calculating f{xj) — ©^^^(j: ) '^i- The 
construction is based on the following observation, which is stated more explicitly 
in [3]. Let v = (wi, ^2, . . . , Wm) be the vector of the function values and let M be 
the mxn incidence matrix of the underlying hypergraph, where the characteristic 
vector of each hyperedge is a row vector of M. If the hypergraph has an empty 
2-core then the linear system M ■ a — v can be solved in linear time. For 
appropriate c this gives expected linear construction time. As before, in the case 
of uniform hypergraphs the space usage is about 1.23 • r bits per key, assuming 
that the values Vj are bit strings of length r. And in our example of mixed 
hypergraphs the space usage is about 1.10 • r bits per key at the cost of a slight 
increase of the evaluation time of /. 

In [8] it is shown how to obtain a retrieval data structure with space usage 
of (1 + e) • r bits per key, for any fixed e > 0, evaluation time 0(log(l/e)), and 
linear expected construction time, while using essentially the same construction 
as above. The central idea is to transfer the problem of solving one large lin- 
ear system into the problem of solving many small linear systems, where each 
system fits into a single memory word and can be solved via precomputed pseu- 
doinverses. As shown in [1] this approach is limited in its practicability but can 
be adapted to build retrieval data structures with 1.10-r bits per key (and fewer) 
for realistic key set sizes. But this modified construction could possibly be out- 
performed by our direct approach of solving one large linear system in expected 
linear time. 

2.3 Perfect Hash Function 

Given a set of keys S, the problem of perfect hashing is to build a function 
h : U — )■ {1,2, ...,n} that is 1-to-l on S. The construction from [3] and [4| 
gives a data structure consisting of a vector a and some mapping that has 
constant evaluation time. Formulated in the context of retrieval, one builds a 
vector V = {vi,V2, ■ ■ ■ ,Vm) such that each key Xj is associated with a value 
f{xj) = Vj that is the index t of the position of a node in the sequence ip{xj). 
This node must have the property that if one applies the peeling process to the 
underlying hypergraph (Algorithm 1) it will be selected and removed because 
it gets degree 1. If c is below the 2-core threshold then with high probability 
for each Xj there exists such an index and the linear system M ■ a = v can 
be solved in linear time. Given the vector a the evaluation of h is done via 
h{xj) = ^{Xj), where l = ®^ev{x,) 
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In the case of a 3-uniforni hypergraph one gets a space usage of about 1.23-2 
bits per key, since there are at most 3 different entries in a. If one appHes a 
simple compression method that stores every 5 consecutive elements from a in 
one byte, one gets a space usage of about 1.23 • 8/5 « 1.97 bits per key. The 
range oi h\& n — 1.23 • m. 

In contrast to the examples above, improving this data structure by simply 
using a mixed hypergraph is not necessarily successful, since the increase of the 
load c is compensated by the increase of the maximum index in the sequence 
ip{xj), which in our example would lead to a space usage of about 1.10-4 bits per 
key for uncompressed a, since we use up to 16 functions for Lp{xj). However, this 
can be circumvented by modifying the construction of the vector v as follows. Let 
G = (5 U {1, 2, ... , n}, E) be a bipartite graph with edge set E = {{a;, gc{xj)} \ 
Xj e S*, t e {1,2,3}}. According to the results on 3-ary cuckoo hashing, see 
e.g. [9,7], it follows that for c < 0.917 (as in our case) the graph G has a left- 
perfect matching with high probability. Given such a matching one stores in v 
for each key Xj the index t of that has the property that {xj, gL{xj)} is a 
matching edge. Now given the solution of M ■ a = v, the function h is evaluated 
via h{xj) = gi{xj) where t = ®i£ip{x )'^i- Since a has at most three different 
entries it follows that the space usage in our mixed hypergraph case is about 
1.10 - 2 bits per key. Using the same compression as before, the space usage 
can be reduced to about 1.10 - 8/5 = 1.76 bits per key. Now the range of h is 
n = 1.10 - m. Solving the linear system can be done in expected linear time. It 
is conjectured that if G has a matching then it is found by the (fc, l)-generalized 
selfless algorithm from [6, Section 5]; this algorithm can be implemented to work 
in expected linear time. 

A more flexible trade-off between space usage and range yields the CHD 
algorithm from [2]. This algorithm allows to gain ranges n — [1 -\- e) ■ m for 
arbitrary e > in combination with a adjustable compression rate that depends 
on some parameter A. For example, using a range of about 1.11 - m, a space usage 
of 1.65 bits per key is achievable, see [2, Fig. 1(b), A = 5]. But since the expected 
construction time of the CHD algorithm is 0{m- (2^ -\- [l/e)^)) [2, Theorem 2], 
our approach could be faster for a comparable space usage and range. 

3 Maximum Thresholds for the Case s = 2 

In this section we state our main theorem that gives a solution for the non-linear 
optimization problem (3) for the case s = 2, that is given two edge sizes we show 
how to compute the optimal (expected) fraction of edges of each size such that 
the threshold of the appearance of a 2-core of a random hypergraph using this 
conflguration is maximal. 

Let fc — (a, b) with a > 3, and b > a. Furthermore, let a ~ {a, 1 — a) and^ 
a G (0, 1], as well as A e (0, -l-cx)). Consider the following threshold function as 
a special case of (2) 

t(A, a, b, a) = ^ . „ . (1 _ e-A)a-i + (l _ «) . 6 - (1 - e-^f-^ ' 



^ We can exclude the case a — 0, since if 3 < a < 6, then it holds that c*(a) > c*(6). 
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We transform i(A, a, 6, a) in a more manageable function using a monotonic and 
bijective domain mapping via z = 1 — e^^ and A = — ln(l — z). Hence the 
transformed threshold function is 



ln(l - z) 
a - a - z"^-! + {1- a) -b - z'' 



T{z, a, b, a) = „ „ , „A . .-1 > (5) 



where z € (0, 1). According to (3) and using T(z,a,b,a) instead of t(A,a, 6, a) 
the optimization problem is defined as 

max min T(z,a,b,a) . (6) 

aG(0,l]2e(0,l) 

For a short formulation of our results we make use of the following three auxiliary 
functions. 

/(,) . -""-f (7) 

g(z.a, b) = J(z} ■ (6 - 1) ■ (a - 1) + — !- + 2 - i. - a (8) 

1 — z 

g . ^g-fe - b - /(^) . (g . (g - 1) . z"-" -&.(&-!)) 

Hz,a,b) = — — — , TT . (9) 

b- {(b-l) ■ f{z) - I) 

Furthermore we need to define some "special" points. 

''^(ir .-/-(^) (10) 

zi — min{z | g{z) = 0} Z2 = max{z | g{z) = 0} . (11) 

It can be shown that if zi and Z2 exist, then it holds z' ^ z\ and z' ^ zi. Now 
we can state our main theorem^. 

Theorems. Let a,b be fixed and letT(z* ,a*) — max min T(z, a). Then the 

^ ^ Qe(o,i];^e(oa) ^ ^ 

following holds: 

1. Let min^ g(z) > 0. 

(i) If h{z') < 1 then the optimal point is (z*, a*) — (z;, 1) and the maximum 
threshold is given by 

T(z*,a*).^Ml^. 



a ■ z 



(a) If h{z') > 1 then the optimal point is the saddle point 



z ,a 



1 



V b/ ' b-a f(z*)-{b-a) 

and the maximum threshold is given by 



T(z*,a*) = -In ( 1 - ''"^ 



For any function (j> — 4>{-,x) we will use (/>(■) and 4>i'yX) synonymously, if x is con- 
sidered to be fixed. 
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2. Let miiiz g{z) < 0. 

(i) If h{z') < 1 then the optimum is the same as in case 
(a) If h{z') G (l,/i(z2)] then the optimum is the same as in case \(ii). 
{in) If h{z') G {h{z2) , h{zi)) then there are two optimal points (z*,a*) and 
{z** , a*). It holds 1/ a* = h{z*) = h{z**) and T{z* , a*) = T{z** , a*). 
The optimal points can be determined numerically using binary search 
for the value a that gives T{zi,a) = T{z2,a), where a is from the 
interval [l//i(zup), l/ft,(zio)] and it holds h{zi) = h{z2) = l/a, with Zi 
from (z;,Zup); o,nd zi from (zioj^r)- (initial) interval for a is: 
• [l/h{zi),l/h{z2)], if zi < z' < Z2, 

. [l/h{z'),l/h{z2)], ifz' < Zl, 
. [l/h{z,),l/h{z')], lfz'>Z2. 

(iv) If h{z') G [h(zi),oo) then the optimum is the same as in case l{ii). 

Sketch of Proof. Assume first that a G (0, 1] is arbitrary but fixed, that is 
we are looking for a global minimum of (5) in ^-direction. Since lim^^o T{z) = 
limz_j.i T{z) = +O0 and T{z) is continuous for z G (0, 1), a global minimum must 
be a point where the first derivative of T{z) is zero, that is a critical point. Let 
z be a critical point of T{z) then it must hold z G [zi, Zr) and a = l/h{z). 

Consider the case mmg{z) > 0. (The case min g{z) = can be handled 
analogously). Since ^^^^ > <^ g{z) > 0, the function h{z) is monotonically 
increasing in G [z;, z^). Furthermore it holds, if g{z) > then z is a local mini- 
mum point of T(z). It follows that for each a there is only one critical point z 
and according to the monotonicity of T(z) this must be a global minimum point. 
Now consider the function of critical points T(z) := T(z, l/h{z)) of T(z, a). It 
holds that 

Vz<z': >0^g(z) >0 and Vz>z': < ^ g(z) > . 

dz az 

It follows that the function of critical points has a global maximum at z' = 

(l) , where z' is at the same time a global minimum of T(z, a) in z-direction. 
If h{z') > 1 then a = l/h{z') G (0,1) and the optimum point (z*,a*) is 
(z', l/h{z')), which is the only saddle point of T{z,a). If h{z') < 1 then because 
of the monotonicity of T(z) the solution for a* is 1 (degenerated solution). Since 
h{zi) = 1 it follows that that (z*, a*) = (z;, 1). 

Consider the case min5(z) < 0. The function g(z) has exactly two roots, 
zi and Z2, and for z G (z/,Zr) the function h{z) is strictly increasing to a local 
maximum at zi, is then strictly decreasing to a local minimum at Z2, and is 
strictly increasing afterwards. Now for fixed a there can be more than one critical 
point and one has to do a case-by-case analysis. A complete proof of the theorem 
is given in Appendix A. □ 

The distinction between case 1 and case 2 of Theorem 3 can be done via solving 
^^^^ = 0, for z G (0, 1), since the function g{z) has only one critical point and 
this point is a global minimum point. Hence, Theorem 3 can be easily transferred 
into an algorithm that determines a* , z* and T{z* , a*) for given k = (a, b). (The 
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pseudocode of such an algorithm is given at the end of Appendix A.) Some results 
for c*{k) = t{\* ,a,b,a*) — T{z* ,a,b,a*) for selected k — {a,b) are given in 
Table 1 and Appendix B. They show that the optimal 2-core threshold of mixed 
hypergraphs can be above the 2-core threshold for 3-uniform hypergraphs. 

4 Experiments 

In this section we consider mixed hypergraphs m a described in Section 1.2. 
For the parameters k = (fci,/c2) G {(3, 4), (3, 8), (3, 16), (3, 21)} and the corre- 
sponding optimal fractions of edge size a* we experimentally approximated the 
point c*(fc) of the phase transition from empty to non-empty 2-core. 

For each fixed tuple {k,a.*) we performed the following experiments. We 
fixed the number of nodes to n = 10^ and considered growing equidistant edge 
densities c = ni/n. The densities covered an interval of size 0.008 with the theo- 
retical 2-core threshold c*{k) in its center. For each quintuple (fci, fc2, a* , n, c) we 
constructed 10^ random hypergraphs „ ^ with nodes {1, 2, . . . , n} and c-a* -n 
edges of size ki and c • (1 — a*) • n edges of size k2. For the random choices of 
each edge we used the pseudo random number generator MT19937 "Mersenne 
Twister" of the GNU Scientific Library [10]. Given a concrete hypergraph we 
applied Algorithm 1 to determine if the 2-core is empty. 

A non-empty 2-core was considered as failure, an empty 2-core was considered 
as success. We measured the failure rate and determined an approximation of 
the 2-core threshold, via fitting the sigmoid function 

ct(c; X, y) ^ {1+ exp(-(c - x)/y)y^ 

to the measured failure rate using the "least squares fit" of gnuplot [21]. The 
resulting fit parameter x = x{k) is our approximation of the theoretical thresh- 
old c*{k). Table 2 compares c*{k) and x{k). The quality of the approximation 
is quantified in terms of the sum of squares of residuals X^ics- '^-'^^ results show 
a difference of theoretical and experimentally estimated threshold of less than 
2-10~^. The corresponding plots of the measured failure rates and the fit function 
are shown in Figures 1, 2, 3 and 4. 
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red data 
(l+e-(--)/^>)-^ 



= 0.00536087 
X - 0.821466 



0-818 0.819 0.82 0.821 0-822 0.823 0,824 0.S25 



red data 



i = 0.00175451 
; a: = 0.851353 



0.848 0.849 0.85 0.851 0.852 0.853 0.854 0.855 



Fig. 1. (fcl,fe2 



(3,4) 



Fig. 2. (fei,fc2) = (3, 
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niea,sureil data 



] Er« = 0.00347648 
i X = 0.910704 



0.907 0-908 0.909 0,91 0.911 0-912 0.913 0,914 



^ lEres = 0.0109112 

J = 0.919848 



0.910 0.917 0.918 0.919 0.92 0.921 0.922 0.923 0.924 



Fig. 3. (fei, fca) = (3, 16) Fig. 4. (fci, fe) = (3, 21) 



(fcl, fc2) 


(3,4) (3,8) (3,16) (3,21) 


C* 
X 

y 


0.82151 0.85138 0.91089 0.92004 
0.82147 0.85135 0.91070 0.91985 
0.00536 0.00175 0.00348 0.01091 



Table 2. Comparison of experimentally approximated and theoretical 2-core thresh- 
olds. The values are rounded to the nearest multiple of 10~^. 



5 Summary and Future Work 

We have shown that the threshold for the appearance of a 2-core in mixed 
hypergraphs can be larger than the 2-core threshold for fc-uniform hypergraphs, 
for each fc > 3. Moreover, for hypergraphs with two given constant edges sizes 
we showed how to determine the optimal (expected) fraction of edges of each 
size, that maximizes the 2-core threshold. The maximum threshold found for 
3 < fci < 6 and ki < k2 < 50 is about 0.92 for fe = (3, 21). We conjecture that 
this is the best possible for two edge sizes. 

Based on the applications of mixed hypergraphs, as for example discussed in 
Section 2, the following question seems natural to ask. Consider the hypergraph 

m a ^^'^ somc fixed upper bound K on the average edge size fc = ^^^j^ • fc^. 

Question. Which pair of vectors k and a. that gives an average edge size below 
K maximizes the threshold for the appearance of a 2-core? That means we are 
looking for the solution of maxmint(A, fe, a.) under the constraint that k < K. 

k,a A>0 
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A Proof of the Main Theorem 

In this section we give the full proof of Theorem 3, i.e. we solve the (transformed) 
non-linear optimization problem (6). As is to be expected, the proof mainly 
employs methods from calculus. 

A.l Preliminaries 

Derivatives. At first we want to determine the partial derivatives of T(z, a, b, a) 
with respect to z and a. To shorten and simplify notation we use the following 
definitions. For all j E N let 

Dj{z,a,b,a) ^a-a-{a-iy ■ z"-^ + {1 - a) ■ b ■ (6-1)^ • z*""^ and 
Zj{z,a,b) = a - {a-iy ■ z"-^ - 6 • (6 - 1)^ • z''^^ . 



The first partial derivatives of T(z, a) are 

dT{z,a) 1 1 In(l-z) Di{z,a) 



dz 1 — z DQ{z,a) z Do{z,ay 

dT{z,a) _ln(l - z) ■ Zo{z) 
da Do{z,ay 

The second partial derivatives of T{z, a) are 

d^T{z,a) _ I 1 2 Di{z,a) 



and (12) 
(13) 

(14) 



{dzf [l-zf Do{z,a) z-{l-z) Da{z,aY 

In(l-z) D2(z,a) - Di{z,a) 2 • ln(l - z) Di{z,a)'^ 



+ 



;2 Doiz,ay z^ Doiz.a)^ 

d^T{z,a) 2 • ln(l - z) • Zo(z)2 



[day Do{z,aY 
d fdT{z,a)\ 1 Zaiz) ln(l - z) Zi{z) 



(15) 
(16) 



dz \ da J 1 — z Dolz.a)^ z £'o(z,Q!)2 

2-ln(l-z) Zo{z) ■ Di{z,a) 
z DQ{z,ay 

Auxiliary Functions. Our analysis is heavily based on three functions, 

/(.) ^ -'""-f (17) 
g{z,a,b)=f{z)-(b-l)-(a-l) + y— + 2-b-a (18) 

1 — Z 

a ■ z-" - b - f{z) ■ (g ■ (g - 1) ■ z"'^ - b ■ {b - 1)) 

b- {(b-i) ■ f{z)-l) 

__ Zo{z,a,b) - f(z) ■ Zi{z,a,b) 
~ fe.z''-i.((6-l)./(z)-l) ' 
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which are shown in Figures 5, 6 and 7 for some parameters a and b. Furthermore 
we make use of the following definitions, 



zi 



1 



a- 1 



1 



6-1 



Zi = min{z | g{z) — 0} 22 = max{z | g{z) — 0} . 



Our line of argument will rely on essential properties of /, g, h and zi, Zr, zi, zi 
and z' . Proving these properties is standard calculus but unfortunately lengthy. 
Therefore the proofs of the next four lemmas are only given in extra sections of 
the appendix. We start with the three auxiliary functions. 

Lemma 1 (Properties of f{z)). 

Let z G (0, 1), then it holds 
(z) fiz) > 1 - z > 0. 
(a) lim^^o/(2^) = 1- 

(iii) lim^^i f{z) = 0. 

(iv) f{z) is strictly decreasing, 
(v) f{z) is concave. 

(w) /(z')>/(^.)-5^. 
{vii) f{z)^-^-2 + b + a. 




Fig. 5. Function f{z) with zi and 
Zr for a — 3,b — 20. 

The proof of Lemma 1 is given in Appendix C. A plot of f{z) is shown in Figure 5. 

Lemma 2 (Properties of g{z)). 

Let S < a < b, z G (0, 1), then it holds 
(i) g{z) is strictly decreasing, reaches a global minimum and is then strictly 

increasing. The global minimum point is the only point where = 0. 

{ii) g{z)>0, Vze (0,z/]. 
{Hi) g{z) > 0, Vz e [zr, 1). 

{iv) Ifimng{z) < then g{z) has exactly two roots, say z\ and Z2, with Zi < Z2 

and Zi, Z2 € (z;, z^). 
{v) Let z > zi then it holds g{z, a, b) > g{z, a,b + 1). 

{vi) For fixed a there is a threshold b' , b' > a + 1, such for a < b < b' it holds 
that min^ g{z,b) > 0, and if b > b' then it holds min^ g{z,b) < 0. 

The proof of Lemma 2 is given in Appendix D. Example plots of g{z,a,b) are 
shown in Figure 6. 

Lemma 3 (Properties of h{z)). 

Let 3 < a < b, z G (0, 1), then it holds 
(i) h{z) has a pole at z — z^. 
{ii) lini2_j.o /i(z) = —00. 
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(iii) Imiz^Zr h{z) — +00. 
{iv) Vz e (0, zi] : h{z) e {-00, 1]. 
{v) Vz € (2:;,-Zr) : h{z) e (l,+oo), and = 1. 
{vi) Vz e {zr, 1) : /i(z) e (-CX), 1). 

^ >0^g(z) >0. 
(■yiw) h{z) is strictly increasing in z £ (0,2;]. 
(ix) h{z) is strictly increasing in z £ [z^, 1). 

(x) If mh\z g{z) > then h{z) is strictly increasing in z £ [zi,Zr). 
(xi) If mhiz g{z) < then h{z) is strictly increasing to a local maximum at zi, 
then strictly decreasing to a local minimum at Z2, then strictly increasing 
afterwards. 

The proof of Lemma 3 is given in Appendix E. Example plots of g{z,a,b) are 
shown in Figure 7. 

60 T n 4 T 

' ~ /i(3.3,10) 
/i(2.3,20) 
/i(2.3,30) 





Fig. 6. Function g{z), Zi,Z2 and 
zi, Zr for a — 3 and b — 20. 



Fig. 7. Function h{z), z\,Z2 
zi , Zr for a = 3 and h = 20. 



and 



Concerning the defined points, we are only interested in how they are related to 
each other. 

Lemma 4. Let 3 < a < b, z G (0, 1), then it holds 

(i) < Zi < Zr < I. 

(a) z; < zi,Z2 < Zr, if zi and Z2 exist. 

(iii) z' € (0, Zr). 

(iv) z' 7^ Zi, z' 7^ Z2, if Zi and Z2 exist. 

The proof Lemma 4 is given in Appendix F. Now we are ready for solving the 
optimization problem. 

A. 2 Analysis 

Assume first that a is arbitrary but fixed, that is we are looking for a global 
minimum of (5) in z-direction. Since 

lim T(z) = lim T(z) = +00 , (20) 

and r(z) is continuous for z € (0, 1), a global minimum must be a point where 
the first derivative of T(z) is zero, that is a critical point. According to (12) 
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critical points in z-direction for unbounded a, i.e., a G M, can be described via 

dTiz) 1 ^ , , -ln(l-z) ^ , , 

dz 1-z z 

^^^^f{z)^a^l/h{z). 

The next lemma identifies and classifies critical points of T{z) for bounded a 
that is for a € (0, 1]. 

Lemma 5. Let a G (0, 1] be arbitrary but fixed. If ^(z) =0 for some z € (0, 1) 
then it holds 

(i) Z € [zi,Zr), 

(ii) if g{z) > then T{z) is a local minimum, 
(in) if g{z) < then T{z) is a local maximum. 



Proof. 

(i) According to (21) we must have a = l/h{z) for a G (0, 1]. Therefore 
it must hold h{z) G (l,+oo). Using Lemma 3{iv) , {v) , (vi) it follows that 

Z e [zi,Zr). 

(ii) Now consider the second derivative of T(z) with respect to z. Ac- 
cording to (14) we have 

d'Tjz) 1 2_ 

(9z)2 (l-z)2 z-{l-z) Do{z) 

ln(l - z) D2{z)-Di{z) 2-ln(l-z) D^jz)^ 

Assume that z G [z;, z^) is a critical point. For the rest of the proof let z = z. 
Utilizing that ~ f(z) it follows that 

d^T{z) 1 2 In(l-z) D2{z) 



{dzY ^°^(l-z)2 z-(l-z)-/(z) ' z2 Dn{z) 
ln(l - z) 2 • ln(l - z) 



z2 • /(z) z2 • /(z)2 



> 



1 /(z) D^iz) 1 

^(1^^ " z-{l-z) ■ D^) + ^ ° 

D2{z) D2[Z) 

Factoring out a from ^j[fy > (1 ^ -z) gives that -Di(z) > (1 — z) • D2{z) is 
equivalent to 

a • (Zi(z) - (1 - z) • Z2(z)) > -& • (6 - 1) • z''~i + (1 - z) • 6 • (6 - 1)2 . z*--! . 
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According to the proof of Lemma it holds a — l/h{z), which can be 
written as 

^ b-z'-^-i{b-l)-f{z)-l) 
Zo{z)-f{z)-Z,{z) ■ 

Division by b ■ z^~^ leads to 

Consider (*). Given that a € (0,1] and z < Zr (Lemma 5(i)) we have, 
according to the definition of and Lemma l(w), that (6 — 1) • f{z) > 1, 
that is the numerator of (★) is larger than 0. Since a > it follows that the 
denominator Zq{z) — f{z) ■ Zi{z) is larger than too. Hence we get 

(az)2 

^{{b - 1) • f{z) - 1) . (Zi(z) - (1 - z) • Z^{z)) > 

((6 - I f • (1 - z) - (6 - 1)) • (Zo(z) - f{z) ■ Z,{z)) 
^Z2{z)-(l-z-{l-z)-{b-l)-f{z)) 

+Z,{z) ■ {{b - If . (1 - z) . /(z) - 1) > Zo(z) . ((1 - z) . (5 - If -{b- 1)) 
^(a-l)2.(l-z-(l-z).(6-l)./(z)) 

+(a-l).((6-l)2.(l-z)./(z)-l)>(l-z).(6-lf -(6-1) . 
Factoring out /(z) • (6 — 1) • (a — 1) gives 

> ^/(z) . (6- 1) . (a - 1) . (6- a) > {b - if - („ - 1)^ _ ^ 

[ozy 1 — z 

^/(z) • (6 - 1) • (a - 1) + + 2- 6- a>0<^ g{z) > . 

(im) Analogous to {ii). 
This finishes the proof of the lemma. □ 



The next lemma can be seen as the central building block for understanding the 
behavior of the threshold function. Using the function g{z) we decide how many 
and which kind of extremal points T(z) has. 

Lemma 6. Let a G (0, 1] be arbitrary but fixed. 

1. Let minz g{z) > then the function T{z) has exactly one critical point z, and 
z € [zi,Zr) is a global minimum point. 

2. Let mhiz g{z) < then there are four pairwise distinct points z^ < zi < Z2 < 
Z2 from the interval [zi,Zr) such that the following holds: 
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(i) For all a with l/a G [1, h{z2)) the function T{z) has exactly one critical 

point z, and z G (z;, z^), is a global minimum point, 
(a) For a with 1/a — h{z2) the function T{z) has exactly two critical points 
zi < Z2, and zi — z^ is a global minimum point, and Z2 = z^ is an 
inflection point. 

(Hi) For all a with 1/a G {h{z2) , h{zi)) the function T{z) has exactly three 

critical points zi < z^ < Z2, and zi, Z2 are local minimum points and Z3 

is a local maximum point, 
(iv) For a with 1/a = h{zi) the function T{z) has exactly two critical points 

zi < Z2, and z\ = z\ is an inflection point, and 2,2 = is an global 

minimum point. 

{v) For all a with 1/a G {h{zi),oo) the function T{z) has exactly one crit- 
ical point z, and z G (z^, z^.), is a global minimum point. 

Figure 8 illustrates the complete case 2 of Lemma 6. The intersection points 
between the function 1/a (horizontal lines) and the function h{z) are the extrema 
of T{z). They are classified depending on the part of h{z) where the intersection 
takes place. 

1.5 



1.4 - 
1.3 - 
> 1.2 



1.1 



1 - 
0.9 -- 

0.7 zi zf 0.75 ^1 0.8 0.85 0.9 Z2 0.95 zl^ 2,. 

z 

Fig. 8. h{z) for a = 3, & = 20, min^ g{z) < 0. 



Proof. 

1. From Lemma 5{i) it follows that all critical points z must be from [zi, Zr). 
Consider the function h{z). According to Lemma 3(f), (x) it holds that for 
each X from [1, +00) there is exactly one z from [zi, z^) such that h{z) — x. 
Furthermore, according to (21) we have ^gf,^-* = a = l/h{z). It 
follows that for each a G (0, 1] there is exactly one z that is a critical 
point, that is it holds a — l/h(z). Let min^ g(z) > then it must hold 
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g{z) > as well. Since z is the only critical point it follows with (20) that 
it must be a global minimum point. 
2. From Lemma 3{xi) we know that for z G [zi,Zr) the function h{z) is 
strictly increasing, reaches a local maximum at zi, is strictly decreasing, 
reaches a local minimum at Z2 and is strictly increasing to +oo afterwards. 
Furthermore it holds g(z) — for z G {21,22} (definition of 21,22), 
g{z) > for 2 < 2i and 2 > 22, as well as 5(2) < for 2 e (21,22) 
(Lemma 3{vii)). 
Consider the condition (21). 
(i) For all a with 1/a G [l,h(z2)) there is, according to Lemma 3, 
exactly one 2 with 1/a — h{z). In addition we have that 2 < 21. 
Utilizing that g{z) > 0, for 2 < 21, (Lemma 2(i), (m)) the claim 
follows by Lemma 5(m). 
(ii) Let 1/a = ^1(22) and let 22 = 22 then according to Lemma 3 there 
is exactly one other point 21, such that a = l/h{zi). Furthermore it 
holds g{zi) > and 3(22) = 0. According to Lemma 3{vii) z\ must 
be a local minimum point. Because of the monotonicity of T(z) (20) 
the other critical point must be an inflection point. Hence 21 is also 
a global minimum point. 
(ill) According to Lemma 3 there are exactly 3 different points 2^, 
1 < z < 3, such that 1/a = h(zi) and ^(^i) = 0, respec- 
tively. Furthermore it holds 21 < 21 < 23 < 22 < 22 and 
g{z\) > 0,(7(22) > 0,(7(23) < 0. From Lemma it fol- 

lows that 2i and 22 are local minimum points of T(z) and 23 is a 
local maximum point of T(2). 
(iv) The case 1/a = l/h{zi) is analogous to the case (ii). 
(v) The case 1/a € {h{zi),oo) is analogous to the case (z). 
This finishes the proof of the lemma. □ 



The last lemma gives a complete characterization of the local extrema of (5) 
in 2-direction including the global minimum for arbitrary but fixed a. It remains 
to find a value a* that maximizes the threshold function at the corresponding 
global minimum in 2-direction. So the point we are looking for could be a saddle 
point of T{z, a). Indeed the following lemma shows that T(z, a) has exactly one 
saddle point for unbounded a, i.e. a € R, and Theorem 3 finally shows under 
which conditions this point is the optimum we are looking for. 

Lemma 7. Let a G M. Then T(z, a) has exactly one saddle point 

(2, a) = (f ) , 51:^ - /(j).(b_a)) ■ 
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Proof. Solving the linear system { ^'^g^'"^ = 0, ^"^g^'"-* — 0} gives 
dT{z,a) 



dz 



= 0^a = l/h{z) 



da Do{z,ay 



<^ a ■ z"-^ = h- z^^^ ^ z = 



There is only one solution of = and according to Lemma 3(z) and 

1 

Lemma l{vi) h{z) is defined at z' = (|) . Hence we get a unique critical 
point {z,a) where z ~ z' and 



6-1 



a = l/h{z) 



h-[lY- ■{J{~z)-{h-l)-\^ 

b-1 
b-a 



b-il)^ ■{.fCz)-{b-l)-l) 

-/(5)-fc'-(f)^-(f-l) 

f{z)-{b-l)-l b-1 



f{z)-{b-a) b-a f{i)-{b-a)- 

To classify this critical point we consider the second partial derivatives 
of T(z,a). We have Zq{z) = and Z\(z) > 0, since 

/flx (a-l)/(fc-a) /a\ 

Z^{z)=a-{a-\)\-^ -b-{b-l)-[-) <0 

It follows that j^£^T{z, a) = as well as g^^ Tl^z, a) > 0. Therefore the 
Hessian matrix H with 

/^(JS^ma) gSsr(5,5)\ /=o >o \ 
[S,ns,a)^Trz,a)) [>O^TCz,a)) 

has determinant det(iJ) < 0, that is (z, a) is a saddle point. □ 



A.3 Putting It All Together 

Now we prove Theorem 3. 
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Proof. Using (21) we can define a function of critical points T{z) of T{z, a) 
as follows 

r(z) :=T(z, 1/hiz)) = b^i 
l/h{z) ■ Zq{z) + b- z" ^ 

_ y^-Zo(z)+ln(l-z)-Zi(z) 

~ 6 • a • (6 - a) • z''+°-2 

The first derivative of T{z) is 

df{z) 1 ( Zo{z) , ln(l-z)-^2(z) 



dz b ■ a ■ {b - a) ■ z''+^-^ \{1 - z)^ z 

b + a-2 ( Zo{z) ln(l - ^) • Zi(z) 



b-a-{b- a) - 
We are interested in the monotonicity of T[z). 
Ofiz) : ^ 



dz 

Zo{z) , ln(l-z)-Z2(z) 
(1-z 



|2 



^ YZ7 + 2) • ^0(2) - • (^2(^) - (6 + a - 2) • Zi(^)) > 

^ _ (6 + a - 2) • Zo(z) + /(z) • Zo{z) • (6 - 1) • (a - 1) > . 

1 — z 

1 

Note that Zo{z) > Q <^ z < (f)''"° = z' . Division by Zq{z) gives, by 
definition of g{z) 

V,>©-:^^<0«,,«>0, 

1. If min^ (/(z) > then according to (★) we have ^'^q^'^ > for all z < z' 

and we have ^g^^-* < for all z > z'. Hence the function of critical points 
has a global maximum in a-direction at z'. Consider the special case 
miriz g{z) = with Zmin = argniin^ ^(z). According to Lemma 2(i) and 
the definition of zi and Z2 we have zi = Z2 = z^i^. From Lemma 4(iw) 
it follows that Zmin 7^ z'- Hence Zmin must be an infiection point of T(z) 
since before and after Zmin the monotonicity is the same. Hence the func- 
tion of critical points has a global maximum in a-direction at z' also in 
this case. 
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(i) If h(z') > 1 then according to Lemma 3(u) we have z' € {zi,Zr). 
It follows from Lemma 6(1) that T(z') is a global minimum in 
z-direction. Hence, (z', l/h(z')) is the optimum point, which is ac- 
cording to Lemma 7 the saddle point, 
(ii) If h{z') < 1 then z' must be from the interval (0, zi] (Lemma 3(iu)) 
and not from the interval (z^,!), (Lemma 3(wi)), since we have 
/(z') > f{zr) (Lemma l{vi)) and /(z) is monotonically decreasing 
(Lemma 1(111)). But if z' < zi then because of the monotonicity of 
T{z) the optimal z value is the nearest feasible critical point. That 
is the optimum point is the (degenerated) solution (zj, 1). 
2. Since min^ g(z) < it follows from Lemma 3(xi) that for z g [z;, z^) the 
function /i(z) is strictly increasing, reaches a maximum at zi, is strictly 
decreasing, reaches a minimum at Z2, and is strictly increasing afterwards. 
Furthermore we have g{z) > for z G [zz,zi), g{z) < for z G (zi,Z2), 
g{z) > for z e (z2, z^.), and (/(z) = for z S {zi, Z2}. An optimal z must 
be global minimum point in z-direction. According to Lemma 5(m) and 
Lemma 6(2) global minimum points are the points from [z;, zi) U (z2, z^). 
(i) An optimal z cannot be from (z2,Zr) since for each z € (z2,Zr) 
there is an e > such that z — e e (z2, z^) and T(z) < T{z — e). 
This descent converges to Z2. But according to Lemma 6(2) Z2 is an 
inflection point and not a global minimum point. Hence the optimal 
z must be from [z;, zi). For each z S (z;, zi) there is an e > such 
that z — £ G [z/, Zl) and T(z) < T(z — e). This descent converges to 

Zl- 

(ii) According to Lemma 4:{iv) z' ^ Z2- It follows that z' e (zj, z^], see 
also Lemma 6(2). An optimal z cannot be from (z2, z^) for the same 
reasons as in case 2(i). 
{at) Consider an arbitrary but fixed a with 1/a e (/i(z2), /i(^i))- Ac- 
cording to Lemma 6(2zii) we have two different points zi,Z2, with 
Zl < Zl < Z2 < Z2, that are local minimum points of the threshold 
function T(z^ a) in z-direction. 

• Let Zl < z' < Z2. Decreasing a (increasing 1/a) by an arbitrary 
small but fixed positive value gives two new local minimum points 
in z-direction, zi H- e, Z2 -I- (5, where e, (5 > 0. According to (*) it 
holds that f'(fi) < f[ix + e) and T{li) > f (is + 5). Hence for 
the left critical point the local minimum in z-direction becomes 
smaller while the potential threshold becomes larger and for the 
right critical point the local minimum in z-direction becomes 
larger while the potential threshold becomes smaller. Increasing 
a by an arbitrary small but fixed positive value reverses the be- 
havior. Assume we have found an optimal a, that is a = a*. 
Decreasing a by some small fixed positive value increases the 
threshold for the left critical point but because of the optimal- 
ity of a we have no global minimum for the left critical point 
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but only a local minimum. Increasing a increases the threshold 
for the right critical point but because of the optimality of a we 
have no global minimum for the right critical point but only a 
local minimum. Hence for a* both critical points z* and z** , with 
1/a* = h{z*) = h{z**), lead to the same minimum in z-direction, 
that is both local minimum points are also global minimum points 
and it holds T{z*,a*) — T{z**,a*) is the optimal threshold. 

• Let z' < z\. Assume that l/a e (/i(z'), ft,(zi)), then a cannot 
be optimal since increasing a by an arbitrary small but fixed 
positive value increases T[z\) as well as T{z2) and one of the 
critical points must be the global minimum point in z-direction. 
Hence the optimum 1/a must be in the interval [h{z2),h{z')]. 

• The case z' > Z2 is analogous to the case z' < zi. 

(iv) According to Lemma 4(iu) z' ^ z\. It follows that z' e [z^, z^). An 
optimal z cannot be from [zj, z\). 

□ 



For given k = {a,b), Algorithm 2 calculates a*,z* and c* = T{z*,a*) of 
Theorem 3 and optimization problem (3), respectively. If one wants to determine 
the optimal values for fixed a but increasing b one can make use of the following 
observation. According to Lemma 2{vi) there is a threshold b' , such that for 
a < 6 < 6' it holds miiiz g{z,b) > and for 6 > 6' it holds g{z,b) < 0. That is 
after reaching 6' we don't need to further calculate the minimum of g{z). The 
following table lists some values for b'. 



a 


345678 9 10 


b' 


16 29 45 62 79 98 117 137 
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Algorithm 2: Optimal Thresholds 
Input: a, b, s (stopping criterion for binary search) 
Purpose: finds optimal thresholds for parameters a and b. 
Prerequisite: suhmutine mmSolveiequation, interval) that returns 

numerical solution of equation within the given interval 

Initialization: 

Zi -h- numSolve(/(z) = z G (0, 1)) 
Zr ^ numSolve(/(2:) = -f^,z € 1)) 

•(-iniinSolve(%^ =Q,z& (0,1)) 
1 

Z' [f] Z\ *r- Z'\ Zi *r- z' 

if g[zg, a, 6) < then 

Z\ ^ numSolve ((/(z, a, 6) = 0, z G (z;, Zg)) 
z-2, *r- numSolve(5(2;, a, 6) = 0, z e (%, -s^-)) 

Oj)timization: 
if h{z\a,b) < 1 then 
I z*^zi; a* ^ 1; T* ^ "'"^^If^ 
else 

if h{z',a,b) < h{z2,a,b) or h{z',a,b) > h{zi,a,b) then 

I ^* ^ a* ^ ^ 7(,.)V^ ; T* ^ In (l - (f ) ^) • (^) 
else 

M ^ zi ; I i— Z2 

if z' < zi then u ^ z' 

if z' > Z2 then / z' 

'^min ^ h{u,a,b) ' "max h{l,a,b) 

while true do 

/>* Qmax— Ctmin I ^ . 

z** -s— numSolve (h{z, a, 6) — ^ = 0, z G (z;, ?i)) 

z* numSolve (h{z, a,b) V = 0, z G {I, z^)) 

t** ^ T{z** ,a,b,a*) 
t* <— T{z* , a, 6, a*) 
if |t* <£then 
I break 
else 

if t* > t** then amin ^ a* 
[_ else amax ^ a* 

return(z*,a*,r*) 
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B Optimal Thresholds 

The following four tables list optimal thresholds for different edge sizes a = ki 
and b = fe, with a e {3, 4, 5, 6} and a < 6 < 50. 
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50 
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6 
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50 


0.98712 4.35225 0.93668 6.91294 0.82405 



Table 3. Optimal values for a = 3 and Table 4. Optimal values for a = 4 and 
a < b < 50. The maximum threshold a ^ b < 50. The maximum threshold 
in this range is about 0.92004. in this range is about 0.82593. 
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93088 


2 


67194 





99807 


6 


04054 





63709 


28 





92783 


2 


62878 





96127 


5 


89081 





70942 


28 





93237 


2 


69378 





99463 


6 


11825 





63717 


29 





92937 


2 


65036 





95979 


5 


96502 





71051 


29 





93379 


2 


71495 





99153 


6 


19490 





63732 


30 





93084 


2 


67130 





95847 


6 


03827 





71163 


30 





93514 


2 


73551 





98873 


6 


27054 





63754 


31 





93223 


2 


69165 





95728 


6 


11061 





71278 


31 





93642 


2 


75549 





98619 


6 


34522 





63781 


32 





93356 


2 


71143 





95622 


6 


18206 





71394 


32 





93765 


2 


77491 





98388 


6 


41900 





63812 


33 





93483 


2 


73069 





95526 


6 


25268 





71512 


33 





93881 


2 


79382 





98178 


6 


49189 





63849 


34 





93604 


2 


74944 





95440 


6 


32250 





71631 


34 





93993 


2 


81224 





97986 


6 


56396 





63889 


35 





93720 


2 


76772 





95361 


6 


39156 





71752 


35 





94100 


2 


83020 





97810 


6 


63523 





63932 


36 





93831 


2 


78556 





95291 


6 


45989 





71873 


36 





94202 


2 


84771 





97648 


6 


70573 





63979 


37 





93937 


2 


80296 





95227 


6 


52751 





71995 


37 





94301 


2 


86481 





97498 


6 


77551 





64028 


38 





94039 


2 


81996 





95168 


6 


59446 





72117 


38 





94395 


2 


88151 





97361 


6 


84458 





64080 


39 





94137 


2 


83657 





95115 


6 


66075 





72240 


39 





94486 


2 


89783 





97233 


6 


91297 





64134 


40 





94232 


2 


85281 





95067 


6 


72643 





72362 


40 





94573 


2 


91379 





97116 


6 


98071 





64190 


41 





94323 


2 


86870 





95024 


6 


79150 





72485 


41 





94657 


2 


92941 





97006 


7 


04783 





64248 


42 





94410 


2 


88425 





94984 


6 


85599 





72608 


42 





94738 


2 


94469 





96905 


7 


11433 





64308 


43 





94495 


2 


89948 





94948 


6 


91991 





72731 


43 





94816 


2 


95966 





96810 


7 


18026 





64369 


44 





94576 


2 


91440 





94915 


6 


98330 





72853 


44 





94892 


2 


97433 





96722 


7 


24562 





64431 


45 





95990 


3 


21627 





94897 


7 


04133 





72973 


45 





94965 


2 


98871 





96640 


7 


31044 





64494 


46 





96531 


3 


36136 





94889 


7 


09535 





73078 


46 





95035 


3 


00281 





96563 


7 


37472 





64558 


47 





96890 


3 


47066 





94885 


7 


14826 





73171 


47 





95103 


3 


01665 





96491 


7 


43850 





64624 


48 





97164 


3 


56291 





94883 


7 


20010 





73252 


48 





95170 


3 


03022 





96424 


7 


50178 





64690 


49 





97386 


3 


64437 





94884 


7 


25089 





73322 


49 





95233 


3 


04355 





96361 


7 


56458 





64756 









































Table 5. Optimal values for a = 5 and Table 6. Optimal values for a — 6 and 
a < b < 50. The maximum threshold a < & < 50. The maximum threshold 
in this range is about 0.73384. in this range is about 0.64823. 
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C Properties of f{z) 

In this section we prove Lemma 1. 

z 

- ln(l - z) > z ^ — ^ — > <4> e"^ > 1 - z / 
I — z 

(a) Applying L'Hopital's rule it follows that 

Um/(.) = lim ^ + . 1 . 

2-)-0 z^O Z z^O 1 

{iii) Applying L'Hopital's rule it follows that 

lim /(z) = hm -^^^^--)-^^--) = lim zMl^ ^ ^ 
= lim 1 — z = . 

df{z) z + ln(l-z) ! , N -, -z y 

V = ^5 ^ < ^ ln(l -z)<-z^l-z<e 

dz z^ 

(„) dV(z) ^ -2z + z2-21n(l-z).(l-z) ! 

(dz)2 (1 - z) • z3 

z^ - 2 • z 
^ — < 2 ln(l - z) 

which is true since it holds 

- Hm^-^o hiz) = lini2_+o f2{z) = and 

__ d/i(z) _ -z^+2-z-2 ^ d/2(z) _ ^-2_ ^ r, 
dz — (l-z)2 ^ dz ^ 1-z ^ 

(vi) First we show that z' is strictly increasing for growing a. Utilizing 
Lemma \ {iv) this implies that /(z') is striclty monotonically decreasing for grow- 
ing a. 

dz^^ , f Ha/b) 1 \ ! 

da ^ \{b-ay^ a-{b-a)) 

b - a f b - a\ b 

<^ > ln(&/a) <:=^ cxp > - 

a \ a / a 



oo 



b-aV 1 



1=0 ^ 



b-a fb-aV 1 fb-a\^ 1 6 , 

a J b a 

>o 
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Now all we need to show is that our assumption holds for the maximum value 
of a, that is a = 6 — 1. 

' b-l\^\^ 1 ^ -ln(l- ^ 1 



6-1 !^ b-1 

b 

-ln(l/6) 1 1 
^^^> — ^ln(5)>l 

which is true since 6 > 4. 

(vii) For z £ (0, 1) we have f{z) e (0, 1). Consider (p{a,b) with 

(p(a, b) = ^ 2 + 6 + a = ^ 1 2 + b + a . 

l-z' 

Let a be fixed. We will show that z' = z'{b) is strictly increasing for growing b. 

dz^^ , f Ha/b) 1 \ ! 

db \ {b-ay b-{b-a)) 

-HI) 1 -h + a (a 

^ (b-a)^^ b-{b~a)^ b ^ U 
4^ X > ln(l + x) , 



for ^ — I + X and — 1 < a; < 0. Hence 
dz' < 

— < ^ X > ln(l + x) ^ X > }^{-iy ■ J— 

i=0 ^ ' 

^x>x -xV2 + x^/i - x^/A / 
^ ' 

<o 

It follows directly that the function Lp{b) is strictly increasing for growing b. 
Consider the minimum of ^p{a + 1, a) = } — l + 2- a = a — 2. Since a > 3 

we have 1^9(0 + 1, a) > 1 which is not in the range of J{z) for z G (0, 1). 

D Properties of g{z, a, b) 

In this section we prove Lemma 2. 

[i) Consider the first derivative of g{z). 

9g(z) ^ ln(l-z)-(b-l).(a-l) ^ 1 ^ (b - 1) ■ (a - 1) ! ^ 

dz (1 — z)^ z 

-/(z)-(b-l)-(a-l) 1 ^ ^ . 

< — • (/(^) - (1 - ^)) 

(6 - 1) • (a - 1) , z , 
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Hence we have RO <^ (b-iy(a~i) ^ ^ ^ —}■ ^ow consider 

gi{z). It holds that 

— lim 51(2) = 0.5, since 

z— i-O 

lim,,(z) = lin.(l^^Hffl^l±^ 

z— >-0 z— s-0 z 

-l.(/(z)-l + z) + (l-z).(^ + l) 

= lim ^ = 1 + lim 

z— !-o 1 z-s-0 dz 

, z + ln(l-z) , 1 + 
= 1 + lim = 1 + Um — 

z^Q z^ z^O 2 • Z 

— 2 . —1 

= 1 + lim = 1 + lim ^ = 0.5 , 

2^o2-z-(l-z) z^o 2 • (1 - z) - 2 • z 

using L'Hopital's rule three times. 

— Iini5i(z) = i=i.(0-1 + 1) = 0. 

z— >1 

— g\[z) is strictly decreasing for growing z S (0, 1), since 

d5i(z) _ -2 • ln(l - z) • (1 - z) + z^ + z^ - 2 • z ^ ^ 
dz z^ 



^2 • ln(l - z) • (1 - z) < z-* + z^ - 2 • z <^ ln(l - z) < 



2.(l-z) 



ff2(z) 



which is true because 
• lim2_i.o ln(l — z) = = lim2_i.o 32 (-z) = and 

, dln(l-z) ^ _^ = -1 - Z < 0. 

dz 1— z dz 

Using that < (^i,_iy(^a-i) follows that for growing z there is a first 

phase with gi{z) > (^i,_iy(^a-i) '^^lidi implies ^^^^<0. Then there is exactly one 
z where gi{z) ~ jj—p^j^^ziYj ' which is a local minimum. After this point we have 
gi{z) < which impilies '^^§f^>0- It follows that the local minimum is 

actual a global minimum. 

(m) If z S (0, Zi] then it holds /(z) > ^j^y > j^y. Furthermore, according 
to Lemma we have > jj^- Let /(z) = and /(z) = as well as 
= with £ > and (5, 7 > 0. Using that /(z) > 0, for z G (0, 1), it follows 

that 

giz) > ^ /(z) . (6 - 1) . (a - 1) + - (a - 1) - (& - 1) > 

^ /(z)2 . (6 - 1) . (a - 1) + (1 + 7) - /(z) • (a - 1) - /(z) • (& - 1) > 
^ (1 + e) • (1 + 5) + (1 + 7) - (1 + e) - (1 + (5) > 
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{Hi) If z e [zr, 1) then it holds f{z) < < Furthermore, accordmg 

to Lemma we have > jj^- Let f{z) = ^^f and f{z) — ^ '^^^^ ^ 
— with e > and (5, 7 > 0. Following the proof of Lemma 2(m) we get 

g{z) > 4^ (1 - e) • (1 - <5) + (1 + 7) - (1 - e) - (1 - (5) > 

(iv) The existance of the roots zi and Z2 follows directly from Lemma 2 (i). 
Moreover, from Lemma 2 (ii), (iii) it follows that if g{z) < for z G (0, 1) then 
it holds z e {zi,Zr). 

{v) Let z > zi, that is f{z) = ^^f for £ > 0. If follows 

g{z, a, b) = fiz) • (fe - 1) • (a - 1) + + 2 - 6 - a 

1 — z 

= f{z) • 6 • (a - 1) + + 2-(fe+l)-a - f{z) • (a - 1) + 1 

= giz, a, + 1) - fiz) • (a - 1) + 1 = g{z, a, 6 + 1) - (1 - e) + 1 
> g{z, a,b + 1) . 

(vi) Assume that there is some b' such that z — miuz g{z, a, b') < 0. Then 
from Lemma 2(ii), (iii) it follows that z > zi. Using Lemma 2{v) we conclude 
that for all 6 > &' it holds that g{z, a,b) < and therefore min^ g{z, a,b) < as 
well. It remains to find one such b'. 

Consider the inequality g{z\ a,b) > which is equivalent to 

(6 - 1) • { fiz')-ia-l) + ^ . -1- -1) - (a _ 1) > . 

gi(o.,b) 

Assume that limt,_>.oo = and limfc_>.oo 32 (fe) < 0. It follows that there must 
be a b' with g{z' {a,b'),a,b') < and thus miiiz g{z' (a, b') , a, b') < 0. 

— limf,_).oo = 0: Assume that it holds limh_^oo -z' = 1- Lemma 1 (iii) gives 
that limfc_^oo f{z'{b)) = = limh^oo gi{b)- 

lim z' = lim exp (in f — ^ ) = lim exp [ — ^— ^ ^ 

6— i-oo 6— >co y V / J fc— >oo y b — a 

/, ln(a) , ln(6)\ / , l/b\ 
= exp lim -j lim = exp — lim = 1 • 

yb— i-oo — a 6— >oo — a / \ 6— i-oo 1 J 

— limb_j.oo g2{b) < 0: Since j^f > iir^i for z G (0, 1), it is sufficient to show that 
limfc^oo f5f7 • 531 = 0. Hence 

lim (2^^0/(^-1), Hm'-^-^ + ^ 



6->-oo 1 — z' fc->oo — 

00 

h^oo U - 1 i 
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Using that %=z'- (^M^ - ^) we get 



b^oo 1-z' 6^00(6-1)2 z' \{b-a)^ b-{b-a) 

= lim — ■( b-1 2 - \ =0- 

6-i.oo z' V (6-a)2 b-{b-a)J 

' . ' 

E Properties of h{z, a, b) 

In this section we prove Lemma 3. 

(i) Consider the denominator of h{z). 

6-((fe-l)-/(z)-l)=0^/(z)= ^ 



1 ' 

which is true for exactly one z from (0, 1), which is per definition z = z^. 

(ii) With lim2_j.o/(2) — 1, Lemma 1(m), and hm2_s.o(l — f{z) ■ (a — 1)) = 
2 — a < —1 we get 

,. a ■ 7^ ■ (1 - f{z) ■{a-l))-b + f{z) -b-jb-l) 

lim h{z) = hm , tt — tt - 00 . 

2^0 z^o 6 • ((6 - 1) • f{z) - 1) 

(zm) It holds f{z) > jj!-j-,Vz G (0, z^). Let /(z) = Consider the limit of 
the numerator of h{z). 

hm a . ^ . (1 - i±f . (a - 1)) - 6 + i±f . 6 . (6 - 1) 

e — 2r 

for some positive constant K (depending on a and b). For the denominator of 
h{z) it holds 

lim 6 • ((6 - 1) • f±f - 1) = 6 • e = +0 . 

Hence lini^_j.^^ ft,(z) = +00. 

(iw) It holds f {z) > ^,Vz e (0,Zi]. Let e > and let f{z) = Hence 

a • z"-'' - 6 - • (a • (a - 1) • z"-'' - 5 • (6 - 1)) ! 

'^'^ = 5-((6-l)-^-l) - ' 

• z"'-'' - (1 + e) • a • z""'' - 6 + (1 + e) • 6 • ^ < -6 + (1 + e) • & • ^ 

• z"-'' - (1 + e) • a • z""'' < 4^ £ > . 

It holds ^ > /(z) > 533;, Vz e {zi,Zr). Let e > 0, and let (5 > with = 
/(z) = i±f . Hence 

a.^-^_6_i=|.a.(a_l).z-^+i±f .5.(5-1)) , 

/i z) = > 1 

5.((fe-l).i±f-l) 

■^a ■ z"-^ - (1 - e) • a • z""'' + (l + 5)-6-6>5-6 

<:^£ • a • z""*" > 4=^ £ > . 
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Note that for z = zi, that is e = 0, we have h{z) = 1. 

(v) It holds f{z) < 5^,V^ e {zr, 1). Let e e (0, 1) and let f{z) = iff. Hence 

a ■ z"-'' - 6 - iff • (a • (a - 1) • z"-'' - 6 • (6 - 1)) i 

"f^' = »-((»-i)-fef-i) " ' 

■ z"-^ • (1 - (1 - £) • ffi) - 6 + 6 • (1 - e) > -e • 6 
- (1 - £) • ffi > <^ £ e (0, 1) . 

We will use the following representation of h{z) 

^ Zojz) - f{z) ■ Z.jz) 

Let hi{z) = (6 - 1) • /(z) - 1. Note that hi{z) ^ 0, if 2; ^ Zr- 
The first derivative of h{z) is 

dh{z) f{z) 



dz b-z^-hiiz) 

(Z,{z) , ^ , (6-l)-(6-l-^) ' 
• (^137 - ^^(^) - (^«(^) - ^(^) • ^f(^^ 

The function h{z) is strictly increasing if and only if 

dh(z 



dz 
dh{z) 



> 

■h,{zf = f{z)-\{h^{z)-{^-Z^{z)) 



dz 

- (Zo(z) - /(z) • Zi(z)) . (6 - 1) . (6 - 1 - ^)] > . 
Note that f{z) is positive for z € (0, 1). So we get 

^ > ^ (4^ - z,{zM{h - 1) • /(.) - 1) > 

This inequality is equivalent to 

^.{h-l)-f{z)-^-Z^{z)-{h-l)-f{z) + Z,{z)> 
1 — Z I — z 



^i(^)..(6-l)./(^) + Zo(^)-(6-l) 



2 



1 - z 



-Zoiz)^^-^-f{z)-Z,{z)-{b-ir 
1 — z 

^f{z) • (6 - 1) • (Zi(^) • (6 - 1) - ^2(2)) + • (Zo • (6 - 1) - 



+ Z2{z)-Zo{z)-{b-lf>Q. 
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Expanding the functions Zj{z) gives 



^ > ^f{z) • (6 - 1) • a • z'^-' ■ {{a - 1) • (6 - 1) - (a - 1)^) 
oz 

+ ^.a.z'^-i.((6-l)-(a-l)) 
+ a-z''-^ -{{a-lf -{b-lf)>0 
^f{z)-{b-l)-{a-l)-{b-a) 



1 



1 - z 



(fe-a) + (a-l)2-(6-l)2 >0 



^/(z) • (5 - 1) • (a - 1) + + 2- 6- a>0 

^g(z) > , 

where we divided by 6 — a which is larger than by definition. 

(vii) This follows directly from Lemma 2(m) and Lemma 3{vii). 
(via) This follows directly from Lemma 2{iii) and Lemma i(vii) 

(ix) For min^ g(z) > this follows directly from Lemma 3{vii). Let min^ g{z) = 
0. According to Lemma 2(i) the point 2;niin — argmin^(7(z) is the only point 
where g{z) — 0. It follows that is the only inflection point of h{z). There- 
fore, according to Lemma 3(wm), h{z) is strictly increasing. 

(x) According to Lemma 2(iv) the function g(z) has exactly two different 
roots zi, Z2 in the interval (z/, Zr) and according to Lemma S{vii) and Lemma 2[i) 
it follows that is strictly increasing for z < zi, strictly decreasing for z with 
Zi < z < Z2, and strictly increasing for z > Z2. Hence the claim follows. 

F Special Points 

In this section we prove Lemma 4. 

(i) follows from Lemma 1(w),(m),(mi). 
(ii) follows from Lemma 2{iv). 
(in) follows from Lemma l{vi). 

(iv) According to the definition of g{z) we have f{z') = g{z') — — 2 + b + a. 
Assume that z' = z\ or z' = Z2, then g(z') — and /(z') = — y^^y — 2 + 6+a 
which is contradiction to Lemma \{vii). 
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